Universal decoherence in solids 
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Symmetry implications for the decoherence of quantum oscillations of a two-state system in a 
solid are studied. When the oscillation frequency is small compared to the Debye frequency, the 
universal lower bound on the decoherence due to the atomic environment is derived in terms of the 
macroscopic parameters of the solid, with no unknown interaction constants. 
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The problem of tunneling and coherence in the pres- 
ence of dissipation has been extensively studied in the 
past, see, e.g., the review of Leggett et al. Q. In this 
Letter we shall revise the case of quantum oscillations of 
a particle (an electron, an atom or a defect) in a double- 
well potential in a solid in the presence of weak dissipa- 
tion due to the non-conducting atomic environment. We 
shall demonstrate that all previous works on this subject 
have missed two important facts. The first fact is that 
the double-well potential formed by the local arrange- 
ment of atoms in a solid is defined in the coordinate frame 
of that local atomic environment, not in the laboratory 
frame. The second fact is that the interactions of the 
tunneling variable with phonons must be invariant with 
respect to global translations and rotations. When these 
facts are taken into account, a simple universal result for 
the decoherence rate can be obtained. 

Let the degenerate minima of the double- well potential 
be located at the points Xq and — Xq inside the solid. In 
the case of weak dissipation, if, at the moment of time 
t = 0, the particle is prepared in the state X — Xq, then 
for any moment t > 



(X(O)X(t)} 



X$e rt cos(uj t) 



(1) 



where (...) means quantum average, Hloq — A is the tun- 
neling splitting, and r <C ojo is the decoherence rate. 
When addressing the effect of dissipation one should an- 
swer two questions: 1) How is A renormalized by the 
atomic environment? and 2) What is the value of T? The 
answer to the first question requires the precise knowl- 
edge of the interactions with the environmental degrees of 
freedom. After the latter are integrated out from the to- 
tal action , A can be computed via the instanton of the 
remaining effective action for the tunneling variable 
For the ohmic interactions, A can be nicely expressed in 
terms of a single measurable dissipation constant p| . In 
this Letter we show that a similar beautiful result follows 
for the decoherence due to the superohmic interactions 
with phonons. In the limit when loq is small compared to 
the Debye frequency, ujd, the decoherence rate L can be 
expressed in terms of measurable constants of the solid, 
with no unknown interaction constants. Similar results 
for problems involving tunneling of the angular momen- 



tum have been obtained earlier |5|, |(j . Here we present 
rigorous derivation of the universal expression for T due 
to phonons for an arbitrary double-well potential in a 
solid. 

The Hamiltonian of the system is 



1 • 2 

U = -mR' +U(Il) + a ik i m R l R k ui m (R) 



( T^P" 2 + KklmUikUlr: 



(2) 



where, R' and R are the radius- vectors of the particle of 
mass m in the laboratory coordinate frame and in the co- 
ordinate frame rigidly coupled with a solid respectively, 
U (R) is the double-well potential, u(r) is the phonon 
displacement field, Uij = ^(diUj +djUi) is the strain ten- 
sor, p is the mass density of the solid, Xikim is the tensor 
of elastic coefficients, and o^zm, ■■■ are coefficients of the 
expansion of the long-wave interaction between R and 
u into a series on Uij . The invariance of the interaction 
with respect to the translation [u(r) — > u(r) + a] and 
rotation [u(r) = $ x r] of a solid as a whole rules out 
combinations of the form Ru and R-Vxu. All interac- 
tion terms, therefore, must be even in the number of the 
spatial components of R. Our result is based, in part, 
upon this observation. Another important observation 
is that the kinetic energy of the particle depends on R', 
while the potential energy depends on R. Substituting 
R' = R + u into Eq. J5J) , one obtains 



n 



U(R) 



-mR 2 

2 

mRu(R) + CtiklmRiRkUlmCR-) 



d 3 r 



[p + m(5(r - R)]ii 2 + \iu m u lk ui m (.3) 



The renormalization of the mass density in Eq. JHJl is to- 
tally insignificant in the long- wave limit (see below), but 
the presence of the effective dynamic interaction, mR-u, 
which is linear in the tunneling variable, is absolutely 
crucial for our argument. 

In this Letter we will treat the renormalized value of 
A as a known parameter that can be obtained from ex- 
periment. We shall start with a symmetric double- well 
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potential. Let the two degenerate minima of U(R) be 
R = ±Ro, so that the ground state and the first excited 
state of the particle are 

|0} = -L(|R ) + |-R )) 

|1) = -L(|R )-|-R )). (4) 

We shall assume that the energy gap, A, between these 
two states, renormalized by the environment, is small 
compared to the distance to other energy levels of the 
particle. Our goal is to compute the decoherence of the 
low-energy states, ) = C% |0) +C2 |1), due to the decay 
of 1 1) into |0), accompanied by the emission of a phonon 
of frequency ujq = A/h. In the limit when ujq <C lojj, the 
phonon wavelength is large compared to the interatomic 
distance, which justifies the use of the long-wave limit 
of the elastic theory. To simplify calculations, we shall 
also assume that the tunneling distance, 2Xq, is small 
compared to the wavelength of phonons of frequency wo- 
Writing R as (X , 0, 0), one has X |±R ) = ±X |±R ), 
while Y I ± R ) = Z I ± R ) = 0. Consequently, X |0) = 
Xq |1) and X |1) = Xq |0), while Y and Z produce zero 
result when acting on |0) and |1). It is now easy to see 
that any combination of the even number of operators X, 
Y, and Z, including X 2 , etc., has a zero matrix element 
between |0) and |1). Thus, the only decohering term in 

Eq. Q is mXiij, with 

(0|1|1) = (0\^(HX - XH)\1) = -iu X , (5) 

where it is explicitly implied that |0) and |1) are the ap- 
proximate eigenstates of the full Hamiltonian, so that 
A = - (0|W|0) is the tunneling splitting renor- 

malized by the interactions. 

We first consider the case of zero temperature. The 
T = decoherence rate is given by the Fermi golden 
rule: 

r = — < «|mRu|j >< j\mR-ii\i > 8{Ei - Ej) 

= ^m 2 \(o\x\i)\ 2 ^\(k,i\l x \o ph )\ 2 s(nu ki -A) (6) 

where \0 p h) and |k, i) are the states of the solid with zero 
phonons and one phonon of the wave vector k and po- 
larization i respectively, and Wki = Cjfc is the phonon 
frequency, with ci being the speed of sound of the polar- 
ization i. The canonical quantization of the phonon field 
yields 

S- = ^= E ^ - «Le-' kr ) ef , (7) 

where V is the volume of the system, <zL and a^i are 
operators of creation and annihilation of phonons, and 



ei are unit polarization vectors. Substitution of Eq. © 
and Eq. into Eq. with account of ((ef) 2 ) = 1/3, 
gives 

1T1R 2 X$U) 2 ^ 

= 3hpV ^ ' ~ ' ' ^ ' 

Replacing the summation over k by V J d 3 k/(2ir) 3 one 
obtains 

6nhp W + cfJ' [> 

where Ct and c/ are the speeds of the transversal and 
longitudinal sound respectively. 

The following observation allows one to understand the 
above result in simple physical terms. Up to a numerical 
factor of order unity, T of Eq. 10 satisfies 

hr~ P 2 j2M, (10) 

with po = mXoujQ being the rms value of the momentum 
of the oscillating particle and M ~ pX 3 being the mass 
of the solid within the volume of dimensions A, where A 
is the wavelength of phonons of frequency wo, averaged 
over polarizations. This volume of the solid adjacent to 
the particle must oscillate together with the particle in 
order to conserve the linear momentum. The latter is 
a consequence of the commonly neglected fact that the 
potential U(R) is formed by atoms that are coupled to 
the rest of the solid, which makes the double-well a part 
of the dissipative environment. The p\j2M width of the 
excited state is, therefore, a consequence of the conserva- 
tion law that mandates the entanglement of the particle 
states with the states of the solid. 

The superohmic case has been studied in Ref. Q, based 
upon the spin-boson Hamiltonian, 

Hsb = -A s x + X s z C a x a 

at 

Here s is the spin-1/2 operator, Ao is the bare splitting, a 
labels coordinates x a , momenta p a , frequencies u> a , and 
masses m a of harmonic oscillators, and C a are constants 
of the linear coupling between the macroscopic variable 
X and the oscillators. The following result was obtained 
for the decoherence rate Q: T = (X 2 /4K)J(A), where 
J{uj) = ^2 a (irC 2 /2m a u! a )5(LL> — ui a ) and A is the renor- 
malized splitting. This, of course, is a correct answer to 
the mathematical problem posed by Eq. I|ll|) . According 
to Ref. Q, for a nonlinear coupling of the form F a (X)x a , 
the quantity C a in the above formulas should be replaced 
by X-^F^Xq) - F a (-±X Q )]. Then, when only even 
powers of X in F a (X) are allowed by symmetry, the re- 
sult for the decoherence rate is zero. As we have seen 
above, in a physical problem of quantum oscillations of a 
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particle or a defect in a double- well potential in a solid the 
linear coupling of the tunneling variable to the phonons 
is prohibited by symmetry, unless one identifies the cou- 
pling with the Galilean transformation term mRii and 
chooses C Q ocoj^. We, therefore, conclude, that in rele- 
vant physical problems the only correct answer for T is 
given by Eq. ©. The beauty of this answer is that it 
does not depend on any unknown interaction constants 
but is expressed in terms of measurable parameters only. 

We shall now generalize our answer to include the case 
of an asymmetric double- well and finite temperature. We 
shall start with the first and will assume that the energy 
minima of U(R) are shifted with respect to one another 
by a small energy e. The origin of the coordinate system 
can always be chosen such that these minima correspond 
to R = ±Ro, with R = (Xq, 0, 0). The relevant two-state 
Hamiltonian of the particle is 



h = — As x + es z . 



(12) 



The solution of the corresponding Schrodinger equation 
yields: 

|0) = i=(CL|Ro> + C + |-Ro» 

|1) = -^(C + |R )-C_|-Ro», (13) 



where 



C± = Vl±e/VA 2 



The energy gap between these states is 



fiuj = VA 2 + e 2 . 



(14) 



(15) 



It is easy to see that the states 111 -it are still the eigen- 
states of X 2 , as for e = 0, so that (0\RiRj\l) = 0, as 
before. Among terms in Eq. which are linear in R, 
the matrix element responsible for the decoherence con- 
tinues to be 



\X\1) 



-icuXoC+C- 



-iuj X . 



(16) 



which is independent of e. The phonon part in Eq. 10, 
however, must be modified by replacing A with the full 
gap, Huj = \J A 2 + e 2 , in the 6- function. At a finite tem- 
perature one should sum up the transitions between the 
initial state, |1), with nki phonons and the final state, 
|0), with nki + 1 phonons, or vice versa for the transition 
|0> -> |1). This gives 



r = 



(17) 



instead of Eq. JHJ. The integration over the phonon 
modes must be accompanied by thermal averaging over 
the number of phonons, with (nki) = [exp^Wki/? 1 ) — 



1] 1 . The final answer for a biased double- well at T ^ 
reads: 



m 2 XQLu 2 uj 3 
6nhp 



2 l\ , [fiuj 
— + — I coth | — 



(18) 



Once again we emphasize that it does not contain any 
unknown interaction constants. The effect of the inter- 
actions has been entirely absorbed into the renormalized 
splitting, Tlusq. 

Among various tunneling problems, the problem stud- 
ied above is relevant to the width of a low-energy opti- 
cal mode that corresponds to quantum oscillations of an 
atom between two isomeric positions in the unit cell of 
a crystal. It describes a kind of an ammonia-molecule 
arrangement of atoms imbedded in a solid. For, e.g., an 
atom of mass m ~ 3x 10 -23 g, oscillating at uj ~ 10 12 s _1 
in a symmetric double- well with X ~ 2xl0~ 8 cm in a 
crystal with p ~ 5g/cm 3 and c t ~ 10 5 cm/s < q at 
T < hu!o, the width given by Eq. I|18() is of the order of 
10 10 s _1 . For tunneling of electrons, the effect is gener- 
ally small due to the smallness of the electron mass. One 
should remember, however, that T of Eq. 1|18|) represents 
the ultimate lower limit of the decoherence rate which is 
mandated by the conservation law. 

To have a complete picture we will show how the above 
treatment of the particle problem can be transformed 
to include the problems of the decoherence of quantum 
oscillations of the angular variable, e.g. spin or an orbital 
moment. This problem is relevant to tunneling of the 
magnetic moment of a molecule between up and down 
directions in a crystal field |5j, |8| • Here we will give its 
general solution in the presence of a bias field. When 
the tunneling variable is the angular momentum L, the 
potential t/(R) in Eq. @ must be replaced with 



7~Lo — PikLiLk + P' ik i m LiLkLiL ri 



(19) 



where j3ik, ^ Hm , etc. are tensors determined by the sym- 
metry of the crystal. The analogy with the particle prob- 
lem is that L in Eq. H19|l is the angular momentum in the 
coordinate frame coupled to the local crystal axes. 

Rotations of L with respect to the crystal axes couple 
to the transversal phonon modes satisfying V-u = 0. 
Consequently, u(R) can be written as JlxR, with 17 
given by f2 = |Vxu. This allows one to write mR-u(R) 
in Eq. © in terms of the angular momentum, 

mRu(R) = (RxmR).O = L O = ^L-(Vxu) . (20) 

While this equation has been derived for L of orbital 
nature, it must equally apply to a spin. One comes to 
this conclusion by considering the effect of the rotation in 
a stationary coordinate frame which axes arc determined 
by the local crystal field at the location of the spin. The 
rotation of that frame due to a transversal phonon is 
equivalent to the magnetic field, which results in the same 
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effective interaction with the spin as it is for the orbital 
moment. 

The other interactions of L with the phonon field are 
of the magnetostriction form, 



(21) 



In the absence of the external field, the symmetry with 
respect to the time reversal requires that all terms in 
Eq. (|21l) contain even number of L-components. The 
external field, H, adds the Zeeman term, — 7LH to the 
Hamiltonian, with 7 being the gyromagnetic ratio. Thus, 
the total Hamiltonian becomes 

H = pikLiLk + p' iklm LiLkLiL m + ... — 7L H 

+ -L-[Vxu(R)] + a ik i m LiLkUi m (Tl) + ... 



F / d\ (ipu 2 + K klM u lm 



(22) 



The main difference of the quantum problem for L from 
the quantum problem for R is that different components 
of the operator L do not commute with each other. For 
certainty, we shall study the situation when the classical 
energy minima of the potential correspond to the classical 
vector L looking in one of the two directions along the 
Z-axis. This could be, e.g., the case of a biaxial crystal 
field and a weak magnetic field applied along the Z-axis, 
H Q = -p z Ll+p x Ll-~fL z H Z) with f3 z ,(3 x > 0. At /3 X = 
the \L) and | — L) eigenstates of L z , would coincide with 
the eigenstates of H.q. Since L z does not commute with 
L x , any small f3 x L x term mixes these states: 



|0) = ±={C.\L) + C + \-L)) 
11} = i(C+|£>-C-|-£» 



(23) 



where C± are given by Eq. JEJ with e = 2-yLH. 

As in the particle problem, the effect of the interac- 
tions in Eq. J22J) on the states |0) and |1) of Eq. lf23|) 
is two-fold. Firstly, interactions renormalize A in the 
expression for the energy gap, Eq. I|15|) . Secondly, they 
generate the finite width of |1) due to the finite probabil- 
ity of the decay |1) — > |0). The latter, as in the particle 
problem, must be considered with account of the sym- 
metry of Eq. (|22|l . In the quantum problem the prod- 
ucts of the L-components in Eq. ll'L'l) must be replaced 
by the symmetric combinations of the Li-operators in 
order to preserve the Hermitian property of the Hamil- 
tonian, e.g., LiLk — ► ^{LiLk + LfeLj), or, equivalently, 
the symmetry of the tensors Pik, etc. with respect to 
the transposition of indices must be enforced. It is then 
easy to see that (0\(LiLk + LkLj)\l) = 0. This is true 
for any combination of the even number of the operators 
Li. Thus, the only term in Eq. i)22[l responsible for the 
decoherence of quantum oscillations between L and — L 
is -|L-[Vxu(R)], which is independent of any interaction 
constants. Consequently, the relevant matrix element of 



the operator of the angular momentum is 

(0|L,|1> = LC+C- = Llu /lu, (24) 

and 

r = -^|( |4|i)| 2 

x |(n ki + l|(Vxu) z |n ki )| 2 <5(w ki - u))) T ,(25) 

where (...)t means thermal average. Substituting here 
the quantized phonon field of Eq. (JJJ), one obtains 



r 



■ coth 



12nhpct 



2T 



(26) 



which is the angular equivalent of Eq. (|18|l . This formula 
generalizes the result of Ref. |5j obtained for H = 0. Same 
as Eq. I|18|) . it contains no unknown interaction constants. 

As in the coordinate tunneling problem, it is interest- 
ing to notice that at H — and T = 0, Eq. I|26l) . can be 
written as TiT ~ L 2 /2I, where / ~ p\ 5 is the moment of 
inertia of the part of the solid of dimensions A~27rc t /cjo, 
adjacent to the particle (molecule) whose angular mo- 
mentum tunnels between L and —L. Thus, the physical 
origin of the decoherence given by Eq. (|26|l is the entan- 
glement of the angular states of the particle (molecule) 
with the angular states of the solid, required by the con- 
servation of the angular momentum. 

In conclusion, we have studied quantum oscillations 
in a double well coupled with a solid. When the oscilla- 
tion frequency is small compared to the Debye frequency, 
the decoherence due to phonons is given by the universal 
formula which contains only directly measurable param- 
eters. It provides the ultimate lower limit for the deco- 
herence rate, mandated by the invariance with respect to 
global translations and rotations. 
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